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Director's  Foreword 


The  US  Army  Research  Laboratory  (ARL)  mission  is  to  “provide  innovative 
science,  technology,  and  analyses  to  enable  full  spectrum  operations.”  As  the 
Army’s  corporate  laboratory,  we  provide  the  technological  underpinnings  critical 
to  providing  capabilities  required  by  our  current  and  future  Soldiers. 

Our  nation  is  projected  to  experience  a  shortage  of  scientists  and  engineers.  ARL 
recognizes  the  criticality  of  intellectual  capital  in  generating  capabilities  for  the 
Army.  As  the  Army’s  corporate  laboratory,  addressing  the  projected  shortfall  is  a 
key  responsibility  for  us.  We  have,  therefore,  identified  the  nation’s  next  generation 
of  scientists  and  engineers  as  a  key  community  of  interest  and  have  generated  a 
robust  educational  outreach  program  to  strengthen  and  support  them.  We  have 
achieved  many  successes  with  this  community.  We  believe  that  the  breadth  and 
depth  of  our  outreach  programs  will  have  a  significant  positive  effect  on  the 
participants,  facilitating  their  journey  toward  becoming  this  Nation’s  next 
generation  of  scientists  and  engineers. 

A  fundamental  component  of  our  outreach  program  is  to  provide  students  research 
experiences  at  ARL.  During  the  summer  of  2016,  we  supported  research 
experiences  at  ARL  for  over  175  undergraduate  and  graduate  students.  Each  of 
these  students  writes  a  paper  describing  the  results  of  the  work  they  performed 
while  at  ARL.  All  of  the  papers  were  of  high  quality,  but  only  a  few  could  be 
selected  for  presentation  at  our  student  symposium.  Several  of  the  presentations  for 
the  selected  research  papers  prepared  this  summer  are  contained  in  this  volume  of 
the  proceedings,  and  they  indicate  that  there  were  many  excellent  research  projects 
with  outstanding  results.  It  is  unfortunate  that  there  was  not  enough  time  for  us  to 
have  all  of  the  papers  presented.  We  would  have  enjoyed  hearing  them  all. 

We  are  very  pleased  to  have  hosted  this  outstanding  group  of  students  for  the 
summer.  It  is  our  hope  that  they  will  continue  their  pursuit  of  technical  degrees  and 
will  someday  assist  us  in  providing  critical  technologies  for  our  Soldiers. 


Philip  Perconti 
Director 
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Introduction 


The  ARL  Summer  Student  Research  Symposium  is  an  ARL  Director’s  Award 
Program  for  all  the  students  participating  in  various  summer  scholarship  and 
contract  activities  across  ARL.  The  goal  of  the  program  is  to  recognize  and 
publicize  exceptional  achievements  made  by  the  students  and  their  mentors  in  the 
support  of  Army  science. 

All  college  undergraduate  and  graduate  students  receiving  research  appointments 
and  conducting  summer  studies  at  ARL  are  automatically  enrolled  in  the 
symposium  program.  As  an  integral  part  of  their  summer  study,  all  students  are 
expected  to  write  a  paper  on  their  work  which  summarizes  their  major  activity  and 
its  end  product. 

The  program  is  conducted  on  two  separate  competitive  levels:  undergraduate  and 
graduate.  The  format  of  the  paper  in  both  levels  is  the  same.  However,  the 
evaluation  will  take  into  consideration  the  difference  in  the  academic  level  of  the 
students. 

All  students  submitted  their  research  paper  for  directorate  review.  Directorate 
judging  panels  selected  one  or  two  papers  from  each  competition  category  for  the 
laboratory-wide  competition  at  the  Summer  Student  Symposium  on  11  August 
2016 

Students  selected  by  their  directorate  for  competition  participated  in  the  one-day 
Summer  Student  Symposium  on  1 1  August  2016.  At  the  symposium,  the  students 
presented  their  papers  to  an  audience  of  ARL  scientists  and  engineers,  including 
the  ARL  Director  and  an  ARL  Fellows  panel. 

This  volume  of  the  Summer  Student  Symposium  Proceedings  contains  many  of  the 
selected  presentations  given  at  the  symposium. 
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Human  Research  &  Engineering  Directorate  (HRED) 
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Applying  Noise  Stimulation  to  Improve 
Trigger  Pull  Steadiness 

Christopher  J.  Blount,  B.S. 

The  University  of  Texas  at  San  Antonio 


What  is  Noise  Stimulation? 

•  Application  of  Gaussian  noise  in  the  form  of 
vibration  to  elicit  an  improved  performance 
response 

•  Imperceptible  (Not  felt  by  participants) 
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How  can  noise  benefit  steadiness? 
Previous  Studies 


•  Past  studies  have  shown  that  introducing  imperceptible 
noise  stimulation  to  biological  systems  can  increase 

force  steadiness.  (Mendez-Balbuena,  2012) 

•  The  input  of  imperceptible  noise  could  decrease  motor 

unit  force  fluctuations.  (Kouzaki,  2012) 
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The  Big  Question: 

How  can  I  determine  an  optimal  imperceptible  noise  level 
to  increase  force  steadiness? 
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Specific  Aim  1 

•  To  define  the  individual  optimal  imperceptible  noise 
stimulation  parameters  that  provide  the  greatest 
improvement  in  force  steadiness  in  the  muscles  that 
produce  the  trigger  pulling  action. 

Specific  Aim  2 

•  To  compare  vibration  stimulus  placements  at  the 
forearm  and  wrist  locations. 


5 


UNCLASSIFIED 


The  Nation’s  Premier  Laboratory  for  Land  Forces 


UNCLASSIFIED 


Methods 

7  male  and  1  female  participants  completed  the 
pilot  study  (n=8) 

•  Participants  were  seated  with  their 
dominant  arm  outstretched  on  a  custom- 
made  apparatus  used  to  simulate  the  grip 
and  trigger  of  an  M4  carbine. 

•  Based  upon  verbal  feedback,  I  determined 
the  imperceptible  vibratory  threshold  per 
individual. 

•  Vibration  levels  were  applied  to  the 
participants  skin  at  the  wrist  and  forearm. 

•  Participants  performed  three  maximal 
voluntary  muscle  contractions  (MVCs). 

•  20%  of  the  greatest  MVC  was  chosen  as 
target  force  during  phase  2,  the  ramp-and- 
hold  task. 
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Ramp-and-Hold  Phase 


Participants  were  instructed  to  trace 
a  line  that  ramps  up  to  20%  MVC 
and  then  maintain  that  force. 

Participants  were  prompted  to  pull 
the  trigger  and  ramp  up  maintaining 
steady  force  for  6  s  and  then  to 
slowly  release  the  trigger. 

Various  imperceptible  noise  levels 
were  applied  to  the  subject  to 
determine  which  amplitude  resulted 
in  greatest  force  steadiness. 
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Data  Reduction  (forearm  placement) 


Root  mean  square  error  (RMSE)  was 
calculated  to  determine  the  error 
between  measured  values  of  applied 
force  and  the  actual  force  trace  during 
the  ramp-and-hold  phase. 

Graphs  were  constructed  using  the 

inverse  RMSE.  (Mendez-Balbuena.  2012) 

Repeated  measures  t-test  was 
performed  to  compare  the  RMSE  with 
zero  noise  (placebo)  and  the  RMSE  with 
optimal  noise. 


Noise  Amplitude  level 


•  For  analysis,  statistical  significance  was 
set  at  0.05. 
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Results  (forearm  placement) 

Paired  t-test  was  performed  to  compare  the  RMSE  with 
zero  noise  and  the  RMSE  with  optimal  noise  values. 


Group 

Placebo 

Optimal  Noise 

Mean 

0.377 

0.239 

SD 

0.1323 

0.0423 

SEM 

0.0467 

0.0149 

N 

8 

8 

(p  =  0.0234) 

Optimal  noise  values  were  significantly  different  than 
Placebo. 
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Analysis  and  Results  ARL 


Wrist  placement 


RMSE  was  calculated  to  determine 
the  error  between  measured  values 
of  applied  force  and  the  actual  force 
trace  during  the  ramp-and-hold 
phase. 

Pearson’s  correlation  was  performed 
on  placebo  amplitudes  to  test  if  a 
learning  effect  had  appeared. 

Repeated  Measures  t-test  was 
performed  to  compare  the  RMSE 
with  zero  noise  and  the  RMSE  with 
optimal  noise. 


Noise  Amplitude  Level 


For  analysis,  statistical  significance 
was  set  at  0.05. 
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Paired  t-test 

•  Paired  t-test  was  performed  to  compare  Optimal  Noise  (ON)  level 
values  from  wrist  placement  and  ON  values  from  forearm  placement 
between  participants. 


Results 


Group 

Wrist  ON 

Forearm  ON 

Mean 

0.273 

0.239 

SD 

0.0436 

0.0383 

SEM 

0.0154 

0.0135 

N 

8 

8 

DF  =  7,  M  =0.018 

p  =  (0.102) 


•  Force  steadiness  in  the  optimal  noise  condition  was  not  significantly 
different  between  wrist  placement  and  forearm  placement. 
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Military 

•  Applying  imperceptible  noise  levels  to  increase  force  stability  could  yield  many 
benefits  to  aid  our  Soldiers  in  tasks  that  require  fine  motor  control. 

•  Further  research  could  be  conducted  on  how  this  technique  can  be  applied  and 
improve  performance  in  other  military  applications. 

•  Shooting 

Civilian 

Using  imperceptible  noise  levels  to  increase  force  steadiness  could  also  benefit 
many  different  areas  in  the  civilian  sector  such  as: 

•  Surgeons 

•  Sports 
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A  Systematic  Study  of  the  Electrical 
Properties  of  Organic  Ballistic  Gelatin 

Ben  Burke 

University  of  Maryland,  College  Park 
Major:  Biological  Sciences 
RSQ-HRED 


•  Background 

•  What  is  electroencephalography?  What  is  a  phantom? 

•  Methods 

•  What  variables  were  tested?  How  were  these  variables  tested? 

•  Results 

•  How  do  these  variables  affect  the  gelatin? 

•  Discussion 

•  What  are  the  next  steps  for  this  research? 
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Electroencephalography  (EEG)  is  the  measurement  of 
brain  activity  via  electric  potentials  at  the  scalp. 


•  EEG  allows  researchers  to  understand  cognitive  effects  of 
various  tasks. 


•  Phantom-electrical  surrogate  with  known  signal  and 
composition  that  can  be  used  to  verify  EEG  devices. 


•  Currently  no  inexpensive,  standard  phantom  recipe  or  material. 


Citation:  https://api.kramesstaywell.com/Content/6066ca30-310a-4170-b001  -  3 

a4ab0 1 3d61  fd/medical-illustrations/l  mages/pnerv201 4033 1 vOQOI _ 
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•  Organic  ballistic  gelatin  can  be  used  to  simulate  human 
tissue  due  to  its  electrical  and  mechanical  properties. 

•  Ballistic  gelatin  is  already  conductive  to  electricity. 

•  No  in-depth  study  of  electrical  properties  of  ballistic  gelatin 

•  The  purpose  was  to  determine  the  extent  to  which  the 
electrical  properties  of  ballistic  gelatin  can  be  controlled. 


Citation:  http://www.thefirearmblog.com/blog/wp-content/uploads/201 5/11/IMGP1 758.jpg 
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Resistance 

Capacitance 

F 

P 

F 

P 

Frequency 

4034.92 

<0.001 

5.441 

0.00482 

Density 

34.93 

<0.001 

1.280 

0.27804 

Frequency*Density 

36.25 

<0.001 

16.468 

<0.001 

A  2-way  ANOVA  was  performed  with  frequency  and  density  as  independent 
variables. 
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Resistance  increased  as  density 
increased 

Noticeably  at  lower  frequency 
levels 


Average  Capacitance  as  a  Function  of  Density 
10 


Density  (percent) 


Inconclusive  effect  of  density  on 
capacitance 

Capacitance  decreased  as 
frequency  increased 
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Key  point:  Density  has  a  relatively  limited  effect. 
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Resistance 

Capacitance 

F 

P 

F 

P 

Frequency 

2402.02 

<0.001 

693.344 

<0.001 

Salinity 

827.79 

<0.001 

560.169 

<0.001 

Density 

11.09 

<0.001 

0.143 

0.7053 

Frequency*Salinity 

451.54 

<0.001 

192.248 

<0.001 

Frequency*Density 

11.92 

<0.001 

1.024 

0.3599 

Salinity*Density 

53.37 

<0.001 

2.993 

0.0307 

Frequency*Salinity*Density 

43.65 

<0.001 

1.099 

0.3621 

A  3-way  ANOVA  was  performed  with  frequency,  salinity,  and  density  as 
independent  variables. 
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Resistance,  10%  Density  Resistance,  20%  Density 


Resistance  tended  to  decrease  as  salinity  increased  across  all  density  levels. 
Capacitance  tended  to  increase  as  salinity  increased  across  all  density  levels. 
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Results  (Warm-Up  Time  Inside  Frame)  ARL 


•  Increase  in  warm-up  time  led  to  inconclusive  effect  on 
resistance  and  increase  in  capacitance. 


Warm-Up  Time  Within  Frame,  Resistance 


Warm-Up  Time  Within  Frame,  Capacitance 
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Results  (Mixing  Time/Warm-Up  Time  Outside 
Frame) 


ARL 


•  Increased  mixing  time  and  warm-up  time  had 
contradictory  effects  on  resistance. 

•  Increased  capacitance  as  mixing  time  and  warm-up  time 
increased. 


Effect  of  Mixing  Time  on  Resistance 

900 


100 

0 

0  s  10  15  20  25  90  JS 

Mlong  time  (minutes) 


Warm-Up  Time  Outside  Frame,  Resistance 
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Effect  of  Mixing  Time  on  Capacitance 


Warm-Up  Time  Outside  Frame,  Capacitance 
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Weighed  a  finished  head  (20%  density,  5%  salinity)  over 
a  16-day  period  to  observe  loss  in  mass 

5.17%  of  original  mass  lost  over  the  period 
Average  loss  of  original  mass  per  day  of  0.3% 


Mass  of  Phantom  Head  Over  Time 
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1  ^  . linear  (Mass  of  Head) 

v  -  -14.422X*  4848.4 

0  2  4  6  8  10  12  14  16  18 

Time  (days) 


The  mass  of  the  phantom 
head  over  a  period  of  almost 
400  h  13 


Density  and  salinity  can  be  changed  to  improve 
resistance  and  capacitance. 

Salinity  has  greater  effect  than  density. 

Increasing  density  to  increase  firmness  can  be 
counteracted  by  increasing  salinity  to  lower  resistance. 

Loss  of  gelatin  mass  over  time  is  an  issue,  but  changes 
to  density  can  be  anticipated. 
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•  Future  research  will  focus  on  utilizing  head  for  testing. 

•  Need  to  look  for  an  alternative,  longer-lasting  material. 

•  Synthetic  gelatin 

•  Carbon  nanofilaments 
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Discussion  (cont.) 


ARL 


Main  takeaways: 

•  Density  and  salinity  can  be  changed  to  tune  electrical  properties. 

•  Organic  gelatin  can  be  used  for  a  phantom  head. 

•  Other  variables  and  materials  need  to  be  tested. 

•  Immediate,  real-world  applications  for  this  research. 

Any  questions? 
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Statistical  Analysis  of  Atmospheric 
Corrections  for  UAS  Acoustic  Signals 


Minas  Benyamin 

University  of  Maryland  at  College  Park 
Electrical  Engineering:  Graduate  2016 


Mentor:  Geoffrey  Goldman 
SEDD/SIP 

Acoustic  and  E&M  Sensing  Branch 
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Unattended  transient  acoustic  MASINT  System  (UTAMS) 
-ARL-developed  system 
-Acoustic  localization  of  weapon  systems 
-Army  top  10  invention  of  the  year:  2004 
-Fielded  in  Iraq  and  Afghanistan 
-Classifier  still  being  developed 


Estimate  target  location 
by  triangulation 
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IT y  RDECOM  New  Classifier  Architecture  ftftL 


UTAMS  signal  processing  Proposed  signal 

architecture  processing  architecture 


Class  I  Gasoline  Powered  UAS  acoustic  data  was  obtained  from 
Webster  Field  using  a  small  tetrahedral  array. 
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Fundamental  frequency  was  calculated  for  each  block  of  time  using 
autocorrelation-based  approach. 

Data  was  resampled  to  have  a  consistent  fundamental  frequency. 


Uncorrected  spectrum 


Corrected  fundamental  frequency 
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RDECOM  Effect  of  Acoustic  Propagation  /ffiL 


Spherical  Spreading 

-  Amplitude  drops  as  a  factor  of  1/Range. 

Frequency  Attenuation 

-  Atmosphere  acts  like  a  low  pass  filter  dependent  on  the 
humidity,  temperature,  and  pressure. 


UAS  Spectra:  Average  Signal  Power  „ 

K  a  a  Range  of  UAS 
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Histograms  of  UAS  for  the 
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Original  Power  PDF 
for  1  to  72.5  m 


Fully  Corrected  PDF  at  Reference  PDF  at 

reference  range  at  24 1  m  rangeof241m 
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Statistical  Test  to  test  for  the  likelihood  that  2  distributions  are  the 


same 


Null  Hypothesis  is  rejected  with  probability  1-a  for  Empirical 
Cumulative  Distribution  Functions  F1  and  F2 


JusteL  A.;  Pena,  D.;  Zamar,  R.  (1997).  "A  multivariate  Kohnogorov-Smimov  test  of  goodness  of  fit". 
Statistics  &  Probability  Letters  35  (3):  251-259. 
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Correction  for  Spherical  Spreading 


Correction  for  Spherical  Spreading 
and  atmospheric  attenuation 


Spherical  Spreading  correction 


Atmospheric  and  Spherical 
Spreading  correction 
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range  intervals 


•  1  to  72.5  m 
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•  Frequencies  of  Interest:  Signal  Harmonics  and  Adjacent  Frequencies  within  ±  5  Hz 

•  Lower  Coefficients  values  indicate  estimated  PDFs  are  more  likely  to  match 

•  Near  Range  <  222  m 
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- .  VT/  Far  Rannp  Corrmarison 


Spherical  Spreading  correction 


•  300  to  485  m 

•  300  to  485  m 

•  485  to  595  m 

•  485  to  595  m 

0.8 

595  to  700  m 

0.8 

595  to  700  m 

!  06 

:  !iii 

l0-6 

W  0.4 
¥ 


0.2 


range  intervals 


•  * 


!iii 

:  .-!  !|'i  :!!:!i!!i! 

Ji  i  *  • !  *  •  i  •  i  •  2 1  5 ! 


;i|||iii|!i!!  ‘  !;!  i!l 

ii'lUHi 


CO  0.4 
¥ 


0.2 


r  : 


Atmospheric  and  Spherical 
Spreading  correction 

range  intervals 


li'i;.  M  In  ,."i 
j!li  !!  ;  • 
;iii  .  i  i ! ! : i H H 1 1 ; | j ; ; ' I ; i » ! ' 

i|  I  l|ip‘!-:  '  i !  ii  iil!S 


1000  2000  3000  4000 

Frequency  (Hz) 


5000 


1000  2000  3000  4000 

Frequency  (Hz) 


5000 


•  Frequencies  of  Interest:  Signal  Harmonics  and  Adjacent  Frequencies  within  5  Hz 

•  Lower  Coefficients  values  indicate  estimated  PDFs  are  more  likely  to  match 

•  Far  Range  >  300  m 
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RDECOM  K-S  Test  Aggregate  Results 


“Average  improvement  in  K-S  value  between  full  correction  and 
only  spherical  spreading 
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psQ  —  refrence  atmospheric  pressure  (atm) 
ps  —  atmospheric  pressure  (atm) 

/  —  acoustic  frequency  (Hz) 

Ta  —  refrence  temperature  (273.16  K) 


T  —  temperature 
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UASs  are  low-cost  asymmetric  threats 

-  Can  easily  be  fitted  with  sensors 
and  even  weapon  systems 

UAS  defeat 

-  Poor  identification  would  lead  to  potential  fratricide 


Classification  of  UASs  remains  a  persistent  challenge 

-  Classifiers  are  often  based  on  feature  matching 

-  Features  can  change  with  propagation  through  the 
atmosphere 
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Joule  Heating  Thin  Film  NiTi 
Cantilevers 

MerricSrour,  2nd  Summer 

Senior  in  Mechanical  Engineering  at  UMBC 
Mentors:  Cory  Knick,  Gabriel  Smith 
August  11,  2016 
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MEMS  Micro-Shutter 

Application: 

•  MEMS  Shutter 

•  Device  Protection 

•  Low  power  on-chip  actuation 

•  Robotics  Actuation 

Passive  MEMS  Shutter  Concept 


Problem:  Micro-Shutterarrav 

•  Requires  small  form  factor,  i.e.  radius  of 
curvature  (~10pm) 

•  Needs  quick  (<100ms)  and  low  power  actuation 

•  Material  properties 

•  Maintaining  shape  memory  effect  in  nano¬ 
scale  films 


Passive 


Active 


My  Focus:  Demonstratingjoule 
heating  active  actuation  method  and 
modeling  device  bending 
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1 .  NiTi  background 

2.  Previous  NiTi  characterization 

3.  Passive  actuation  methods 

4.  Active  joule  heating  actuation  method  (my  summer 
research) 

5.  Modeling  radius  of  curvature 

6.  Summary  and  future  work 


Nickel  Titanium  (NiTi)  is  a  Shape  Memory  Alloy  (SMA),  meaning  it  is  a  type  of 
metal  alloy  that  can  recover  seemingly  permanent  strain  when  heated. 

SMAs  and,  specifically,  NiTi  have  2  phases  -  martensite  and  austenite. 

•  Martensite  :  Low  temperature  phase  (ductile) 

•  Austenite  :  High  temperature  phase 


Martensite 


Austenite 


Tools  to  characterize  NiTi  : 

•  Differential  Scanning  Calorimetry  (DSC) 

•  X-ray  Diffraction  (XRD) 

•  Wafer  Bow  Tool 

•  Stress  vs.  Temperature  using 

E  h. 


Stoney's  equation:  a  = 


\R  R0J 


6(1— v)  h 

Kahn  H,  Huff  MA,  Heuer  AH.  The  TiNi  shape-memory  alloy  and  its  applications  for  MEMS.  J.  Micromech.  Microeng.  1998;8:213-221. 
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Previous  characterization  of  NiTi  material 


XRD  Data  Comparing  Varying 
Thicknesses  of  NiTi 


l  M 


—1.6  Jim  NiTi 
-527  nm  NiTi 
-341  nm  NiTi 
-270  nm  NiTi 


DSC  Full  Heat/Cool  Cycle 


•  Primary  martensite  peak  :  41.5° 

•  Primary  austenite  peak  :  42.5° 


Reversible  phase 
change  above  RT  -58°C 


Previous  characterization  of  NiTi  material  : 


Stress  vs.  Temperature  for  NiTi  on  Si 


3000 


20  40  60  80  100 

Temperature  (°C) 


•  Completely  reversible 

•  Large  recovery  stress 

•  Low  residual  stress 


Reversible 

Large  recovery  stress 
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Device  Fabrication  &  Operation  ARL 


When  released,  cantilevers  curl 
up  due  to  stress  in  bimorph 
layer  (Pt  /  NiTi). 

•  CTE  mismatch,  and  cooling 
from  600°C  deposition  or 
anneal  temperature 

When  NiTi  is  actuated,  it  goes 
from  martensite  to  austenite 
phase,  doubling  in  rigidity,  and 
goes  to  its  shape  remembered 
state. 


1-  Pt  evaporation  and  NiTi  sputter 

2-  NiTi  wet  etch,  and  Pt  ion  mill  using  5214  pattern 

3  -  Xenon  Difluoride  dry  release  from  Si  wafer 

4  -  Thermal  actuation  (laser,  joule  heating,  etc.) 
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Laser  Actuation  Demonstration 


0.000  s 


Laser  Intensity  vs. 
Actuation  Time 
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Previously  had  actuated  the  cantilevers  using  heater  and  laser, 
now  actuating  the  devices  using  joule  heating. 


i 


Device  Design  in 
AutoCAD 

•  3  arrays  with  different  width 
cantilevers 

•  20|jm  ,  15|jm  ,  10|jm 


Advantages  of  joule  heating  : 


•  Increases  potential  on  chip  applications 

•  MEMS  Sensing /actuating 

•  Switching 

•  Low  power 

•  Ability  to  resonate  actuators 

•  Find  thin  film  NiTi  limits 
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NiTi  Wafer  Recipes 
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Device  Recipes 


NiTi 

Thickness 

(nm) 

Deposition 

Temperature 

(°C) 

Post 

Anneal 

(600°) 

600 

600 

No 

300 

600 

Yes 

300 

600 

No 

300 

400 

Yes 

300 

400 

No 

*AII  NiTi  was  sputtered  on  a  20  nm 
layer  of  Pt 


Reason  for  trying  300  nm 
@  400  °C 


•  Reversible  300  nm  NiTi 
film  on  Si 
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Microcombustion  and 
Thermophotovoltaics 


Noah  Zecher-Freeman 
Power  Components  Branch 
Mentor:  Mike  Waits 
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Soldier  Power 


ARL 


Soldier  is  overburdened  with  weight,  restricts  mission  length. 
Current  lithium  batteries  have  a  specific  energy  of  160  Wh/kg. 
For72-h  mission,  20  W  Soldier  power  requires  9  kg  of  batteries. 
Thermophotovoltaic  (TPV)  cells  take  advantage  of  the  specific 
energy  of  hydrocarbons  to  reduce  the  weight  carried  by  Soldiers. 


selective  emitter 


Microreactor 


electrical  output 


7  —  ?!  thermal'  i^l  spectral'  Vvf'  ^IPVcell) 

=  50%  •  (70%  ■  25%)  =  8.8%  FY16  BPP  Objective 
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RDECOM  Previous  Microreactors  MKL 


DSI  reactor  used  to  study  the  role  of 
confinement 


Tolmachoff,  et  al.,  Applied  Energy,  2014 


Frontiers  MR  2.0  (FY16) 

Adds  heat  recirculation 
Validate  reaction  &  CFD  models 
Operate  in  vacuum  chamber 


Frontiers  MR  1.0  and  1.1  (FY14-FY15) 

Studied  performance  with  high  heat  loss  &  short  residence  times 
Visualize  gas  reactions 

Develop/validate  reaction  models:  surface  +  gas-phase  fuel  decomposition 
Tolmachoff  ef  a/  .  Combust.  Flame,  2015 


Frontiers  MR  2.0  flow 


Previous  modeling  detailed  how  to  maximize  efficiency  of  reactor 


Jirr-r- 
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Radiation  to 
PV  cell 


Energy  balance: 

Combustion 
Desired  surfaces 
Heat  lost  to  exhaust 
Undesired  surfaces  Catalyst  Coating  / 
Conduction  along  channels 


Radiation  to 
rTPV  cell 


4.23*02  5  73*02  7.22*02  8.72*02  1.02*03  1.17  *03  1.32*03  1.42*03 


[1]  Tolmachoff  et  al..  Combustion  and  Flame,  162  (2015),  pp  3674-3680. 

[2]  ED.  Tolmachoff  et  al.,  Proc.  Western  States  Section  of  the  Combustion  Institute  Fall  Meeting.  2015. 
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Schematic 
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pGC 

Air  and  Fuel 
Flow  Controllers 


heating  element 


Thermocouple  and  Pyrometer  readings 


20  mm  x  20  mm  Reactor  Surface  Temperature  (e=.4) 
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Reactor  Length  (mm) 

20  mm  x  20  mm  Reactor  Surface  Temperature  (e=.9) 


Hthermal  53% 

Peak  Temperature  =  1405K 
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15 

1000-1050 
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Reactor  Length  (mm) 


Hthermal  53% 

Peak  temperature  =  1148K 
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Channel  1.  lOin 
MS5A  Heated 
Injector.  Backflush 
Results 

Pk  *  Same 


Retention  rune 


Area 


Concent/ arion 


4  Oxygen 
8  Nitrogen 

-2  i.fohanf 


1 10  Carbon  Monoxide 


0.529 

0.714 

mu 


1628 


17290 

24113046 

:iuii 


600648 


0.000 
732136.072 
925  629, 


20594.298 1 


Totals 


24752324 


“’53635.999 


C  hannel  2.  10m 
PPU  Heated 
Injector.  Backflush 
Results 


Pk* 

Same 

Retention  Time 

Area 

Concent /  arion 

h 

Carbon  Dioxide 

0509 

5789016 

147M4.77S  | 

3 

ttnyiene 

U.535 

7556 

20  '.323 

4 

Ethane 

0552 

11766 

289.732 

Acetylene 

0.000  BDL 

6 

Propane  Propylene 

0805 

8151 

87.004 

7 

Unknown 

1.131 

178882 

0.000 

Totals 

36194938 

906037.610 
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60% 
|  58% 

OJ 

is  56% 

LU 

|  54% 
|  52% 
50% 


C12H26  +  18.5 02  ->  12C02  +  13 H20 


(p  =  1 ,  stoichiometric 
(p  <  1,  lean  mixture 
(p  >  1,  rich  mixture 


60  80  100  120  140  160  180 

%  Theoretical  Air 

Figure  5.  Combustion  Gas  Concentrations  at  Percent 
of  the  Theoretical  Combustion  Air 


TSI  Incorporated,  Combustion  Analysis  Basics 
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Carbon  Mass  Balance 

100.00% 

95.00% 

90.00% 

CM 

8  85.00% 

80.00% 

75.00% 

70.00% 

0.8  0.85  0.9  0.95  1  1.05  1.1  1.15  1.2  1.25  1.3 

Equivalence  Ratio 

— *—  Run  1  —♦—Run  2  — *—  Run  3 


•  Thermal  efficiency  (radiated  heat  from  desired  surfaces) 

•  Over  5  runs,  thermal  efficiency  was  51  -58% 

•  Demonstrated  20%  TPV  efficiency  [1]  combined  with  50% 
combustion  efficiency  for  total  efficiency  of  10%  allowing  for 
energy  density  of  ~700  Wh/kg,  2.1  kg  for  72  h,  20  W  mission 

•  Combustion  efficiency 

•  Carbon  balance  suggests  that  combustion  is  nearly  complete 

•  Heat  lost  to  exhaust 

•  ~3.5%  energy  lost  is  due  to  high-temperature  exhaust 

•  Radiated  heat  from  undesired  surfaces 

•  >15%  energy  lost  is  due  to  radiated  heat  from  undesired 
surfaces,  which  can  potentially  be  reflected  back 

•  Heat  lost  to  conduction  along  channels 

[1]  Wemsman  et  at,  IEEE  Transactions  on  Electron  Devices,  51  (2004),  pp.  512-515. 
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Future  Efforts 
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Emissivity  of  the  reactor  is  unknown  (with  higher  emissivity  while 
oxidized,  decreases  when  reduced  inside  vacuum);  the  diffuse 
reflection  complicates  measurements. 

Current  reactor  has  a  bulge  in  it  due  to  being  heated  and  high 
pressure  inside  the  reactor  and  low  pressure  outside;  earlier 
modeling  done  is  now  inaccurate. 

New  reactors  coming  in  are  to  be  reinforced. 

Experiment  with  complex  fuels  rather  than  simple  fuels. 

Determine  improved  methods  of  depositing  catalyst;  transition  to 
other  catalysts  for  better  resistance  to  fouling  from  sulfur. 


Bill  Allmon 
John  Little 


UNCLASSIFIED//For  Official  Use  Only 


The  Nation’s  Premier  Laboratory  for  Land  Forces 


Approved  for  public  release;  distribution  is  unlimited. 


52 


M 


UNCLASSIFIED//For  Official  Use  Only 


US.  ARMY 

FSDECOM 


ARL 


Back  up 


Tolmachoff  et  al.,  Combustion  and  Flame,  162  (2015),  pp.  3674-3680. 


V=  3  SLPM  (156  W),  Twall  =  1300  K 
surface  (solid  lines)  gas  phase  (dashed  lines) 


UNCLASSIFIED//For  Official  Use  Only 


The  Nation’s  Premier  Laboratory  for  Land  Forces 


Approved  for  public  release;  distribution  is  unlimited. 


53 


XL 


UNCLASSIFIED//For  Official  Use  Only 


US.  ARMY 

RDECOM 


Selected  Results:  Effect  of  Flow  Ra  *ARL 


Tolmachoff  et  al.,  Combustion  and  Flame,  162  (2015),  pp.  3674-3680. 


Tvvall  =  1300  K,  26=0.5  mm 
2, 3,  4  SLPM  =  104, 156, 208  W,  respectively 


/  1* 
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Visualizing  Ferret  Brain  Neural  Pathways  to 
Understand  Mechanisms  of  Blast  Brain  Injury 


Skylar  Bodt,  CQL  Intern 

Mentors:  Autumn  Kulaga  &  Karin  Rafaels,  ARL-SLAD 


RDRL-SLB-WJuly  28th,  201 6 


Fallston  High  School  (2015) 

•  US  Congressional  Art  Competition 

•  GEMSatARL  2012  &  2013 


University  of  Maryland,  College  Park  (2019) 

•  Major:  Cell  Biology  and  Molecular  Genetics 


•  Activities:  Beta  Psi  Omega  Professional  Biology  Fraternity,  UMCP 
Scholar,  Quality  Enhancement  Systems  and  Teams  (QUEST)  Honors 
Program 


Future  Career  Aspirations 

•  Graduate  degree  in  Medical  Illustration 

•  Continue  research  as  a  biomedical  illustrator 
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Project  Goal 


ARL 


GOAL:  Further  TBI  research  through  visual  reconstruction  of 
neuroanatomy 


Objectives: 

1.  Hand  render  general  images  of  a  ferret  brain  in  ventral,  lateral,  and  dorsal  views,  clearly 
identifying  regions  of  the  brain.  Modify  iteratively  using  Adobe  Photoshop  and  Illustrator. 

2.  Within  these  images,  locate  the  position  of  whole  cerebrum  sections  taken  from  a  ferret 
brain  and  stained  with  B-amyloid  precursor  protein  (B-APP)  or  a  neurofilament  stain  (NF), 
both  intended  to  show  axonal  damage. 

3.  Illustrate  neural  pathways  based  on  the  2D  stain  images. 
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Current  threats,  such  as  Improvised  Explosive 
Devices  (lEDs),  make  Soldier  exposure  to  blast  a 
continuing  concern. 

Advances  in  soldier  protective  equipment  have 
increased  soldier  survival  rate,  but  created  an 
increase  in  nonlethal  Traumatic  Brain  Injury  (TBI). 


This  survival  rate  increase  brings  into  focus 
nonlethal  TBI,  and  a  need  to  understand  the 
mechanism  of  blast  injury  to  the  head  and  neck 
region  (1). 


Studying  soldier  TBI  risk  due  to  blast  effects  is 
vital  for  both  the  treatment  of  injury  and 
development  of  protective  equipment  (1). 
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Why  is  it  Important? 

•  The  ferret  brain  is  important  in  research  because  it  helps  bridge  the  gap 
between  rat  and  human  brain  research. 

•  It  has  more  similar  surface  area  to  a  human  but  maintains  the  size  of  a 
small  rodent  brain.  (1) 

fgkl 

What  are  the  challenges? 

•  There  is  little  research  available  on  the  ferret  brain, 
so  there  are  gaps  in  knowledge. 

•  Hardly  any  visualizations  have  been  created,  so 
as  a  result,  known  information  is  difficult  for  researchers  to  interpret. 


Important  Previous  Research 

•  An  experimental  program  in  2010  used  a  ferret  animal  model  to  explore 
brain  damage  after  exposure  to  blast  overpressures  in  a  controlled 
environment.  (1) 

•  I  used  this  research  as  the  foundation  for  my  own  project. 


UNCLASSIFIED 


The  Nation’s  Premier  Laboratory  for  Land  Forces 


Approved  for  public  release;  distribution  is  unlimited. 


59 


UNCLASSIFIED 


Approved  for  public  release;  distribution  is  unlimited. 


60 


UNCLASSIFIED 


★  US.  ARMY 

IT j  RDECOM  Methods:  Objective  One  ARL 


UNCLASSIFIED 


The  Nation’s  Premier  Laboratory  for  Land  Forces 


Approved  for  public  release;  distribution  is  unlimited. 

61 


UNCLASSIFIED 


★  •  US  ARMY  _  ^  - 

IT y  RDECOM  Methods:  Objective  One  ARL 


UNCLASSIFIED 


The  Nation’s  Premier  Laboratory  for  Land  Forces 


Approved  for  public  release;  distribution  is  unlimited. 


62 


UNCLASSIFIED 


U.S.  ARMY 

RDECOM 


Methods:  Objective  One  ARL 
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Dorsal  View  ARL 


Olfactory  Bulb 


Longitudinal  Fissure 


Orbital  Gyrus 


Cerebellum 


Posterior  Ectosylvian  Gyrus 


Dorsal  view 


Posterior  Sigmoid  Gyrus 


Coronal  Gyrus 


Anterior  Sigmoid  Gyrus 


Lateral  Gyrus 


Anterior  Ectosylvian  Gyrus 


Suprasylvian  Sulcus 
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Ventral  view 


|  Olfactory  Bulb  Optic  Nerve  Mamillary  Body  *>ons  j|||  Medulla 

Lateral  Olfactory  Tract  Infundibulum  Crus  Cerebri  |H|  Pyramid  Spinal  Cord 
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Lateral  view 


|  Olfactory  Bulb 

■ 

Coronal  Gyrus 

■ 

Lateral  Gyrus 

■ 

Posterior  Sigmoid  Gyrus 

■ 

Anterior  Ectosylvian  Gyrus 

■ 

Posterior  Ectosylvian  Gyrus 

■ 

Suprasylvian  Gyrus 
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•  Determined  Bregma  (0-value)  on 
rendered  images  (2) 

•  Bregma:  a  location  on  the  front 
of  the  skull  where  plates  come 
together 

•  Difficult  to  determine  in  ferrets 

•  Deciphered  section  placement 
from  description  of  #  mm  from 
Bregma  on  the  ferret  brain  sample 

•  Adjusted  scale  to  reflect  -  3.2  cm 
ferret  brain  depicted  on  10.4  cm 
rendered  image,  (-325%  increase) 

•  Uploaded  image  into  Adobe 
Illustrator  to  show  slice  placement 


Ferret  skull  and  brain 


|iumiii|uiimu  tniuHt|iumiut|  mm 

-20  -10  Bregma  10  20 


A 


•  36  sections  coronally  cut  over  whole  cerebrum,  20  B-App,  16  NF  (1) 

•  Vague  description  where  slices/cuts  were  made  and  slide 
numbering  from  anterior  to  posteriori  NF  &  1 -B-APP  /  mm(1) 

•  Slide  placement  estimate  by  Talairach  Stereotactic  Coordinates 


NF  JU  B-APP 
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il  i-S.miY 


Dorsal  view 

Zero  marker  for  slices 
(Bregma) 

Positive  slice  distance  -* - |— ►  Negative  slice  distance 
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All  Images  Created 
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Lateral  view 


Dorsal  view 


Ventral  view 


Lateral  view 


Dorsal  view 


Ventral  view 
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Notional  View  of  Neural  Pathways 
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•  Understanding  TBI  is  of  critical  importance  to  developing  better 
protective  equipment  and  smarter  assessment  of  that  equipment. 

•  Visualization  of  the  ferret  brain  and  regions  will  help  researchers 
identify  important  structure  related  to  TBI. 

•  Determining  slide  section  locations  for  NF  and  B-APP  stains 
allows  a  visual  reference  that  links  the  neural  pathway  damage 
seen  on  the  sections  to  a  specific  region  within  the  brain. 

•  Completed  work  in  depicting  connections  of  neural  pathways 
across  sections  will  move  toward  a  3D  understanding  of  TBI. 

•  Future  work  might  improve  the  visualization  to  allow,  for  example, 
“removable”  section  images  from  the  main  image  or  3D  rendering. 


•  I  appreciated  this  opportunity  to  be  a  Medical  Illustrator  Intern 
withARL/SLAD. 

•  The  biology,  mainly  anatomy,  required  forced  me  to  expand  my 
knowledge. 

•  I  now  know  much  more  about  the  regions  of  the  brain  and  the 
coordinate  system  used  to  locate  structures. 

•  I  had  not  previously  done  anything  at  this  level  with  Adobe 
Photoshop  and  Illustrator;  I  had  only  experimented.  This  summer 
forced  me  to  become  much  more  proficient  in  both,  but  I  have 
much  to  learn. 

•  I  want  to  thank  my  mentors  Autumn  Kulaga  and  Karin  Rafaels  for 
guiding  me  this  summer  in  Illustration  and  Biology. 
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Stains 
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B-amyloid: 

•  Beta  amyloid  is  a  protein  fragment  detached  from  an  amyloid 
precursor  protein  (APP). 

•  Increased  B-APP  expression  is  associated  with  brain  injury  (6). 
Neurofilament: 

•  Neurofilaments  are  important  elements  in  the  cytoplasm  that 
supports  the  cytoplasm  in  the  axon. 

•  Staining  may  give  insight  to  axonal  damage  in  the  brain  (6). 
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Early-Stage  Damage  in  CFRP  using 
Active  Thermography 


Student:  Micheal  Wallace 
Mentor:  Dr.  Mulugeta  Haile 
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•  Motivation 

•  Theory 

•  Experimental  Setup 

•  Results  and  Discussion 

•  Summary 

•  Future  work 
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•  The  overarching  objective  is  to  discover  in-situ  sensing  technique 
for  detection  of  early-stage  fatigue  damage  or  “damage 
precursors”  in  polymer  composites  using  the  materials  local 
thermal  properties. 


Damage  Precursors 
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Micro-cracks 

formation 


Local  heating 


Dislocation  pile-up 
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current 
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0.8 
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Macro-crack 

formation 


Cracks 

delam/debond 

Change  in  dynamics 
Change  in  structural  response 
Severely  degraded  module  efficiency 
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Service  Life 
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•  Fatigue  damage  results  in  microscopic  changes  in  materials 
(dislocation,  crazing  and  slip  bands). 

•  Microscopic  changes  in  turn  results  in  changes  in  the  behaviorof 
lattice  vibration  (or  Phonon). 

•  In  nonmetals,  heat  conductivity  is  mainly  due  to  phonons. 

•  All  other  things  constant,  changes  in  materials  thermal  property, 
such  as  conductivity,  diffusivity,  heat  capacity  may  indicate  early 
stage  fatigue  damage. 


Early  stage 
material 
degradation 
due  to 
Fatigue 


Microscopic 

changes 

such 

crazing, 

micro- 

craking 


Changes  in 
lattice 
vibration 
(Phonon) 


Thermal 

property 

changes 
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Data 

acquisition 


Flash  Source  IR  radiation 


240v,  20A 


Specimen 


Flash  Source 
240v,  20A 


ID  Heat 
conduction 


Steps: 

1 .  Heat  the  specimen  surface  with 
high  intensity  flash. 

2.  Measure  surface  temperature 
change  AT  as  a  function  of 
time,  i.e.  AT  =  f(t). 

3.  Process  the  temperature  time 
(T-t)  data  to  detect  abnormalities. 


Principles  of  active  thermography 
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Transient  Heat  Conduction-  ID 


Flash  Source  IR  radiation 


240v,  20A 


Specimen 


Flash  Source 
240v,  20A 


ID  Heat 
conduction 


Principles  of  active  thermography 


Governing  equation  (Non- 
dimensional  form 

—=v'-e 

8t 


d(X, 0)  =  L  ^  =  0  at X=0  and  r>0 

-=-Bi.e 

sx 

Solution 

0{x’t)^=Aie'K'cos{^f-j  T>02 

Where 

6{X  t)  =  ^  [X.  t)  -  r.  Dimensionless  temperature 
X 

X=  —  Dimensionless  distance 
L 

Dimensionless  heat  transfer 
Bi  =  k  =  —  coefficient.  Biot  Number 
k 

at 

t  =  —  Dimensionless  time,  Fourier  Number 
L‘  - 
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Army 

Platform 


Principal  component  - 
Fatigue  sensitive  panel 


Geometry 

Max 

Min 

R-ratio 

Spectrum 


W=  25mm,  t=  3.5mm,  A=  86  mm2,L=  167  mm 
Pm*  =  22.2  kN,  Om*  =  260  Mpa 
Pm*  =  2.2  kN,  Om*  =  26.0  Mpa 
0.1 

Uniform 


Constant  spectrum  load 

800 

700 

600 

I500 

I400 

£ 

W  300 
200 

100 

0 

06  08  1  1.2  14  16  18  2 

Strain  (m/m)  xio  3 


Ultimate  Tensile  Strength  *  720  Mpa 

— 


36%  UTS  -260  Mpa 
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Experimental  Set-Up  ARL 


MTS  with  Specimen 


IR  Camera  with 
EcoTherm  Hood 
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RDECOM  Experimental  Procedure  MRL 


Specimen 
are  subjected 
to  cyclic 
fatigue 
loading 


Active 

thermographic 
images  are 
then  captured 
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IR  Images  at  0  Cycles  IR  Images  at  80K  Cycles 
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Cycles 


g> 
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- 0  Cycles 

- 100  Cycles 

- 200  Cycles 

- 300  Cycles 

- 400  Cycles 

- 500  Cycles 

- 750  Cycles 

-  1000  Cycles 

-  1250  Cycles 

1500  Cycles 
2000  Cycles 
—  2500  Cycles 
3000  Cycles 
3500  Cycles 
4000  Cycles 
4500  Cycles 
5000  Cycles 
6000  Cycles 
7000  Cycles 
8000  Cycles 
9000  Cycles 
10000  Cycles 
12500  Cycles 
15000  Cycles 
17500  Cycles 
20000  Cycles 

- 25000  Cycles 

- 30000  Cycles 

- 35000  Cycles 

- 40000  Cycles 

- 45000  Cycles 

- 50000  Cycles 

- 60000  Cycles 

- 70000  Cycles 

- 80000  Cycles 
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Interpreting  the  Log  Plot 


y-intercept 
indicates  the 
amount  of  energy 
absorbed. 


In  (time) 


Simple  Interface:  Adiabatic  Boundary 


steel 

|  1  cm 

heat  transfer  through  back  wall  (adiabatic) 


The  time  t*  characterizes  the 
disruption  of  the  diffusion 
process  by  an  internal  interface. 

In  the  adiabatic  case,  t*  is  the 


intersection  of  the  straight  lines. 
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Results 
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Second  Derivative:  A  Sensitive  Indicator  of  V 


d2ln(T-T;) 
d2  ln(t) 


i\ 


ln(t*) 


ln(t) 


t*  signals  the  arrival  of  the 
heat  pulse  at  the  back  wall 
interface 


Derivatives  are  highly 
sensitive  to  t*  transition 
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I^T)  RDECOM  Discussion 

ll  U.S.ARMY 


10 


4 _ i _ i _ i _ i i i i— i— i _ i _ i _ i i i i (—i—i _ i _ i _ i i i i  i  r 

2  3  4  5 

10  10  10  10 

Fatigue  Cycles 


Fatigue  Cycles 
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Aerodynamic  Optimization  of  a 
Supersonic  Bending  Body  Projectile  by  a 
Vector  Evaluated  Genetic  Algorithm 


Justin  Paul 

Mentor:  Dr.  Sidra  Silton 
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About  Me  AftL 

•  I  was  born  and  raised  in  Austin  Texas  with  a  strong 
interest  in  mathematics  and  physics. 

•  I  am  currently  undergraduate  senior  in  Aerospace 
Engineering  at  The  University  of  Texas  at  Austin. 


® TEXAS 

The  University  of  Texas  at  Austin 


Interests  include  rocketry,  machine  learning,  and 
juggling 
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I.  Internship  Objectives 

II.  Project  Description: 

a)  Research  purpose 

b)  Projectile  studied 

c)  Aerodynamics  considered 

d)  Conditions  considered 

III.  Implementing  Optimization 

IV.  Simulation  Methods  Explored 

V.  CFD  Methodology 

VI.  Genetic  Algorithms: 

a)  Description 

b)  Advantages  and  Disadvantages 

c)  Methodology 

VII.  Current  Results: 

a)  Analysis 

b)  Possible  conclusions 

VIII.  Path  Forward 

IX.  Improving  Projectile  Shape: 

a)  Discuss  more  realistic  projectile  shapie 
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1.  Internship  Objectives 

ARL 

Primary: 


•  Study  ballistic  flight  dynamics  and  bending  body  literature 

•  Learn  to  use  and  experiment  with: 

-  Missile  DATCOM 

-  Solidworks  Flow  Simulation 

-  CFD++ 

-  MIME 

•  Optimize  control  parameters  for  a 
bending  body  projectile 


of  Grav  ity) 


Secondary: 

•  Study  evolutionary  programming 

•  Write  genetic  algorithm  to  complete  optimization  (CML-VEGA) 
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II.  Research  Purpose 


ARL 


Maneuvering  flight  bodies  for  small  diameter 

munitions  research  area 

-munitions  that  are  both  accurate  and  precise 

-minimize  collateral  damage 

Current  control  authority  is  limited: 

-moving  control  surfaces 
-thrust  vectoring 

Current  research  done  on  deflectable  projectile 
structure: 

-focused  on  deflecting  the  nose  cone  at  the  base  of  the 
nose  cone 

Novel  aspects  of  this  investigation: 

-multiple  bends 

-bends  located  on  the  body 


Projectile:  Air  Force  Firmer 


-Landers  et  al.  (2003)  -Shoesmith  et  al.  (2006) 
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II.  Aerodynamics  Considered 


ARL 


Wind  Vector  (M^) 


-  AOA 


£ 


M 


^n 


•  Pitching  Moment  Coefficient:  Primary  characteristic  of  control  authority 

M 

C»  -  I - 2 - 

sd 

•  Axial  Force  Coefficient:  Characteristic  of  adverse  aerodynamics 

c  -  Fx 
2  KPooMoo  s 

•  Normal  Force  Coefficient:  Secondary  characteristic  of  control  authority 

Fn 


Cu  — 


1  2 
2  7 Poo  Moo  s 
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AOA  =  Angle  between  the  trajectory  vector  of  the  projectile  and  the 
freestream  vector 

Mo,  =  The  freestream  Mach  number 
Poo  =  Ambient  pressure  in  the  freestream 
Tm  =  Ambient  temperature  in  the  freestream 


AOA  =  {-8: 1:8}° 

AOA  (deg) 

MACH  2 

MACH  3 

MACH  4 

-8 

Case  1 

Case  18 

Case  35 

Mm=  {2:1:4} 

-7 

-6 

Poo  =  101325  kPa 

-5 

-4 

Tm  =  293.15  K 

-3 

-2 

51  Total  Cases 

-1 

Case  8 

Case  25 

Case  43 
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III.  Implementing  Optimization  ARL 


Initialization 

(pseudo-random 

population) 


Geometry 
(defined  by 
control 
parameters) 


Legend 

MIME  &  CFD++ 

SOLIDWORKS 

CML-VEGA 


GD 


Convergence 

check 


Generate 
aerodynamic 
coefficients 
(through  CFD) 


Individual 
Optimization 
(Cm*,  CN*,  CX*) 


Optimized 

Individual 
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Aero-prediction  code 


Computational  Fluid  Dynamics  Solvers 


Missile  DATCOM? 

-  Fast  analysis 

-  Semi-empirical 
(limited  accuracy) 

-  Limited  modeling 
capabilities 


Legend 

Advantages 

Disadvantages 


SOLIDWORKS  Flow 

Simulation: 

-  Integrated  into  CAD 
software 

-  Fast  meshing 

-  Limited  licensing  and 
resources 

-  Limited  to  windows 
systems 


CFD++: 

-  More  sophisticated 
solver 

-  Near  unlimited 
licensing  and 
resources 

-  Compatible  with  UNIX 
workstations  (and 
therefore  HPC) 

-  Requires  MIME 
software  for  meshing 
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V.  CFD  Methodology 

ARL 

Version  15.1.1  ofCFD++: 

•  Steady-state  and  quasi-steady  sweep 
RANS  simulations 

•  Realizable  kc  turbulence  model 

•  Finite  volume  framework 


Version  5.1  of  MIME: 

•  Unstructured  tetrahedrals  with 
prism  boundary  layer 

•  Local  refinement  on  nose  cone  tip,  bends  and  edges 

•  Curvature  refinement  across  nose  cone 

•  Global  mesh  size  typically  9.1  to  9.6  million  ceils  (taking  advantage  of 
longitudinal  symmetry) 
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Defining  Aspects  of  a  Genetic  Algorithm: 

•  Metaheuristic  (practical  solutions) 

•  Stochastic  (search  and  propagation  utilize  pseudo  random  number 
generation  techniques) 

•  Conditionally  Global  (initial  population  size,  crossover  probability  and 
mutation  rate  dependency) 

Advantages: 

•  Robust  composition  is  insensitive  to  initial  condition  variation 

•  Relatively  simple  to  implement  on  relatively  complex  systems 

Ex:  size(a1)*size(  f>1)*size(4>2)  =  16*16*16  =  4096  configurations 

Disadvantages: 

•  Does  not  scale  well  -  increasing  quantity  of  control  parameters  increases 
the  diversity  of  design  space  exponentially 

•  High  dimensional,  multimodal  problems  require  sophisticated  and 
expensive  fitness  functions 
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VI.  Genetic  Algorithms 
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The  Model 


Optimize  for  maximum 
surface  area 


•  ©  • 

Fitness  Evaluation 

•  • 

•  • 

and 

•  • 

Natural  Selection 

Initial  Population 


No 


Fit  Population 


Convergence? 


Yes  ^ 


Recombination 


Offspring  Population 


Optimized  Individual 


Selected  Parents 


•  • 


•  • 
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For  a  given  set  of  aerodynamic  coefficients:  min:  niche :  max 

max  —  min 

size  =  - - - 

niche 

component  score  =  priority  coefficient  *  value  score 

N 

final  fitness  score  =  y  (component  score) j 
7=i 


value 

value  score  =  -  *  size 

max 


<(1/2)Cm,(1/3)Cx,  (1/6 )C„> 

( final  fitness  score)i  Wj 

fitness  proportionate  probability  =  - 7 - — - = - 

maximum  final  fitness  score  wmax 


Minimization  modification:  value  score  =  — 7 —  *  size 

value 

*The  complete  objective  function  includes  CFD  simulation  calculations 
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VII.  CM  vs  AOA  ARL 

Pitching  Moment  Coefficient:  Primary  characteristic  of  control  authority 


3 


AOA(°) 
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Mach  2  ,  AOA  8° 


Generation  1 


Generation  2 


Individual  4 


Individual  1 


Individual  26 


Individual  4 


Individual  29 


Individual  6 


UNCLASSIFIED  -  APPROVED  FOR  PUBLIC  RELEASE 

15  The  Nation’s  Premier  Laboratory  for  Land  Forces 

UNCLASSIFIED 

FkWT)  RDECOM 

{!■«»— »]* 

VII.  Generation  Statistics 

ARL 

length(0:0.5:7. 5)  =  lengthfl. 25:4. 5:68.75)  =  16 


Mach  =  2 


Generation  1 

Mean  Std.  Dev.  4.71 


Mean 
Std.  Dev. 


Ql  ^1 

8.97  8.19  8.10 

4.71  4.89  4.52 


Generation  2 


Mean  Std.  Dev.  3.69 


22%  reduced  standard  deviation 


Angle  of 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

AOA 

Attack 

Mean 

-8° 

-7° 

-6° 

-5° 

-4° 

-3° 

-2° 

-1° 

0° 

1° 

2° 

3° 

4° 

5° 

6° 

7° 

8° 

Standard 

Deviation 

3.77 

3.95 

3.85 

4.21 

2.82 

2.99 

2.88 

2.55 

4.63 

4.57 

4.34 

4.54 

3.43 

2.21 

4.13 

4.27 

3.60 
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•  Algorithm  functionality  is  not  completely  confirmed  -  observed  increase 
in  performance  will  verify  optimization. 

•  Central  Objective:  complete  optimization  of  unique  problems 
(characterized  by  CM  vs  AOA) 

Ex:  AOA  =  8°  for  Mach  =  {2:1:4} 

•  Complete  statistical  analysis  of  algorithm  convergence  to  understand  why 
some  problems  may  converge  faster  than  others. 

•  Look  for  significant  differences  between  problems  with  related  conditions 
to  derive  conclusions. 

•  Analyze  flow  field  characteristics  of  optimized  individuals  to  understand 
increased  performance  compared  to  the  baseline  straight  body 
configuration. 
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IX.  Improving  Projectile  Shape 

ARL 

Use  Bezier  control  points  to  create  a  continuous  shape  (Li  et  al.  2016) 


Fourth  Order  Bezier  Curve 
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•  Natural  Selection: 

-median  population  fitness  is  calculated 

-individuals  at  or  below  the  median  fitness  are  removed  from  the  population 
to  create  the  “fit”  population 


•  Parental  Selection: 


1 

N 

1 

N 

1 

N 

1 

N 

uniform  probability 

■ 

w± 

wmcoc 

w2 

wmcuc 

w3 

wmax 

Wi 

wmax 

* 

fitness  proportionate  probability 

selected  parent 

if - —  >  x,  0  <  x  <  1 

Wmax 

individual  i  is  the  selected  parent 
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Three  dimensional  cut-plane  of  fin  region  cells 
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Comparison  and  Optimization  of  Solid- 
State  Joining  Techniques  for  Dissimilar 

Metals 


Miriam  Silton 

Coatings,  Corrosion,  and  Engineered  Polymers  Branch  (CCEPB) 
Mentor:  Dr.  Robert  Jensen 
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□  About  me 

□  Background 

□  Goals 

■  Fabricate  and  characterize  standardized  joints 

■  Produce  baseline  performance  predictions 

□  Conclusions  and  Future  Work 
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Friction  Stir  Welding  (FSW) 


AD  (Accumulating 
direction) 


Close-up  of  welding  area 

Ultrasonic  Additive  Manufacturing  (UAM) 
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Adhesive 


Machine  lap 
joints  and 
thread  butt 
joints 


Clean  and 
pretreat 
bonding 
surface 


Apply  adhesive 


Cure  with  heat 
and  applied 
pressure 


Sand  fillet 


Cold  Spray 


Thread,  insert 
steel  pin,  and 
bevel 


Clean  bonding 
surface 


Spray  while 
rotating 


Wait  until  cool 


Machine  off 
extra  spray 
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>  Standard  tensile  testing  of  butt  joints  (ISO  6922/ASTM  2095) 

■  Adhesive:  dissimilar  Al  alloys 

■  Cold  Spray:  dissimilar  Al  and  Mg  alloys  with  Al  spray 

>  Post-corrosion  tensile  testing  of  adhesive  lap  joints 

■  30  days  in  salt  fog  chamber 

■  5%  NaCI-water solution  @  35°C  (ASTM  B117) 


>  Cyclic  fatigue  testing  of  adhesive  lap  joints  (ASTM  D3166) 
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>  Cold  spray:  high  potential  vs.  large  standard  deviation 

>  Adhesive:  consistent  +  comparable  strength 


80  AI2024/6061  + 


0  0.2  0.4  0.6  0.8  1  1.2  1.4 

Extension  (mm) 
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>  Current  Literature 

■  Hybrid  joints  combining  adhesives  with  UAM  and  FSW 

■  Cold  spray  is  new  technique  with  limited  research 

>  Hybrid  Cold  Spray-Adhesive  Butt  Joint 

■  High  strength  of  cold  spray 

■  Reliability  of  adhesive 

■  Optimal  jointfor  this  geometry? 

•  Hypothesis:  Combination  will  eliminate  crack  and  reduce  stress 
concentrators  at  edge. 

>  Considerations 

■  Ensure  weld  heat  doesn’t  degrade  adhesive 

•  Degradation  temp  ~360°C 

•  Surface  during  cold  spray  application  ~60°C 
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“Note:  Hybrid  tensile  test  was  paused  for  about  5  seconds  around  100  MPa  and  data  was 
spliced  together  to  create  data  curve  shown. 
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Hybrid  Analysis 
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Cohesive  failure 
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>  Decreased  strength  and  displacement  retention  with  adhesively 
bonded  dissimilarjoints  -  what  happens  with  hybrid? 


20 


Corrosion 


Extension  (mm) 


30  day  salt  fog  exposure 
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>  Will  hybrid  joint  improve  fatigue  resistance? 
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Cyclic  fatigue 
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>  Conclusions 

■  Geometrical  limitations  of  solid-state  techniques 

■  Corrosion  has  noticeable  effect  on  dissimilar  joints 

■  Opportunity  for  hybrid  joints  to  increase  mechanical  performance 

>  Short-term  Future  Work 

■  Fabricate  and  test  hybrid  butt  joints 

■  Fabricate  hybrid  lap  joints 

■  Develop  FSW  and  UAM  methods 

>  Long-term  Future  Work 

■  Fatigue  and  corrosion  analysis 

■  Full-scale  testing 

■  Modeling 
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QUESTIONS??? 
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Tool 
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Transducer 


* 

Frictic 


AD  (Accumulating 
direction) 


Metal  tape, 


VD  (Vibrating 
direction) 


(Rolling  direction) 


Close-up  of  welding  area 


Powder  Feeder  Supersonic  Nozzle 
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>  80/20  bar  with  spring  screw  that  can  be  tightened  to  apply  pressure 
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>  .005”  diameter  glass  beads  to  ensure  even  bond  line  thickness 

>  Chopped  fibers  for  increased  viscosity  and  strength 


•8& 
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Shear-Induced  Initiation  of  Nitromethane 


1Brad  A.  Steele  and  2James  P.  Larentzos 
1Ivan  I.  Oleynik 

University  of  South  Florida  Physics  Department 
Tampa,  FL 

2Army  Research  Laboratory,  WMRD  Aberdeen 
Proving  Ground,  MD 

DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited 


The  Nation's  Premier  Laboratory  for  Land  Forces 


UNCLASSIFIED 


Energetic  Materials  ARL 


Trinitrotoluene  (TNT) 


Reaction:  A - ►  C 


A  better  understanding  of  the  initiation  of  energetic  materials  can 
make  them  safer  and  allows  for  more  advanced  energy  release. 

The  initiation  is  critically  dependent  on  the  reaction  pathway  and  the 
energy  barrier  height  and  width. 
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Mechanical  Initiation  ARL 


•  Shock  compression: 


V  =  Vn 


Energetic 

Material 


Hydrostatic  compression: 


Energetic 

Material 


Uniaxial  compression: 


Energetic 

Material 


Force  F 


Shear  stress: 


Energetic 

Material 


Area  A 


Shear 

stress 

F 

T~A 
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Nitromethane 


ARL 


Nitromethane:  (N02)(CH3) 

■  A  prototypical  EM 

■  Widely  studied  and 
experimentally  well- 
characterized 

■  A  model  for  understanding 
chemistry  in  larger 
explosives  with  similar 
functional  groups 

■  No  phase  transitions 
reported  in  hydrostatic 
pressure  experiments 


Nitromethane 
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Shear-Induced  Initiation  ARL 


Reacts  at  just 

18  GPa 

under  shear! 

Reacts  at  30 
GPa  under 
compression 
alone 


250 


750  1250 

Wavenumber  (cm1) 


1750 


Experimental  collaborators  at  ARL  (Tim  Jenkins)  compress  and  shear  NM 
in  a  rotational  diamond  anvil  cell  (RDAC)  and  measure  Raman  spectrum. 


The  Nation's  Premier  Laboratory  for  Land  Forces 


Computationally  investigate  the  chemical 
reaction  mechanisms  that  occur  during  the 
shear  stress  and  initiation  of  nitromethane. 


■  Develop  a  theoretical  methodology  that  can  accurately 
describe  the  response  of  nitromethane  under  shear 
stress. 

■  Investigate  the  role  crystallographic  orientation  plays 
on  the  shear  sensitivity  of  NM. 

■  The  initiation  of  NM  may  be  related  to  the  presence  of 
grain  boundaries  in  the  crystal  -►  investigate  the  shear 
dynamics  with  grain  boundaries  in  the  crystal. 
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•  Atoms  follow  Newton's 
laws  via  classical  molecular 
dynamics. 

•  A  REACTIVE  potential  is 
used  to  model  atomic 
interactions. 

•  Potential  shows  good 
agreement  to  experimental 
EOS  and  melting  lines. 


•  Crystal  orientations  and 
grain  boundaries  built 
using  code  developed  at 
ARL. 
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Fully  periodic  crystals  built  with  various  orientations 
sheared  along  all  three  directions 

■  Initially  relaxed  at  20  GPa  and  298  K,  then  sheared  (shear  rate  =  -0.386°/ps) 


<01-1> 
orientation 
sheared 
along  xy 
angle 


90° 


W'wm  --  - 

f  UNCLASSIFIED 


shear  direction 


xyz.=90° 


xy  z.  =65.77° 
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Significant  orientation  dependence 


•  Shear  stress  along  sheared  direction  shows  a  gradual  decrease  to  -2  to 
-5  GPa,  then  a  sharp  increase,  accompanied  by  a  transition  to 
disordered  NM  state. 


•  Size  of  shear  stress  and  shear  angle  are  highly  dependent  on  the 
orientation. 


<01 -l>  Oriented  along  z-axis 


—  xy  angle  deform 

—  xz  angle  deform 

—  yz  angle  deform 


70  80 

Angle  (degrees) 


90 


Increasing  shear 


Disordered  state 


61.22° 


xyz.  =61.22° 


<0 1  - 1  > 


out  of  plane 
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<01-1>  Orientation  sheared  along  xy^i 

■  When  sheared  along  this  direction,  NM  molecules  can 
rotate  slightly  into  empty  space  perpendicular  to  shear 
direction. 
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<01-1>  Orientation  sheared  along  yzz 

■  NM  molecules  once  again  rotate  to  be  aligned  along  y- 
axis. 

■  Plane  of  NM  molecules  shift  slightly  along  shear 
direction  (y-axis). 


•  Several  hundred  bi-crystal  grain  boundaries  (GBs)  are 
created. 

•  Atom  positions  and  crystal  lattice  relaxed  at  20  GPa 
and  300  K. 

•  GBs  are  ranked  according  to  GB  energy: 

g  _  U GB  ~~  U perfect  *  ^ atoms 

2  *  Area 


•  Structures  with  lowest  GB  energy  are  chosen  on  which 
to  perform  shear  simulations. 
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Grain  Boundary: 
<011>/<011> 
along  z 

<0  27  -14>/< 0 
27  14>  along  y 


u 


Shear  Grain  Boundary  AKL 


Angle  90 


SB^»V2»2<V2’2/»  f*5’Sf* 
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Notice  that: 


Grain 

Boundary 


Grain 

Boundary 


1.  The  entire  grain  rotates  with  the 
shear. 

2.  The  molecules  do  not  rotate  until 
a  very  large  shear  angle  is 
reached. 

3.  Disordered  NM  originates  at  GB 
and  grows. 
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QQ  COO 


Tintra  Avg.  (K) 


Tintra  Avg.  (K) 

450 


78.25° 


Syz=-0.58  GPa 
T=289  K 
P=19.86  GPa 


Syz=-1.32  GPa 
T=340  K 
P=19.93  GPa 


67.80° 


Tintra 


(K) 


58.70° 


Tintra 


(K) 


Syz=-1.63  GPa 
T=378  K 
P=20.09  GPa 


Syz=-1.62  GPa 
T=424 K 
P=19.55  GPa 
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Shear  Along  xz 


ARL 


Grain  Boundary: 
<011>/<011> 
along  z 

<0  27  -14>/<  0  - 
27  14>  along  y 
*Shear  direction 
into  plane 


Angle  90 
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Grain 

"Boundary 


Grain 

Boundary 


Grain 

Boundary 


Notice  that: 


1.  Transformation  to  disordered 
state  occurs  suddenly 
throughout  the  entire  material. 

2.  Stress  still  grows  because 
transformation  is  incomplete  at 
~  80°. 

3.  Much  different  than  yz 
direction. 


F^WT) 
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Future  Work 


ARL 


Chemical  reactions  have  not  occurred 
during  shear  simulations! 

•  Perform  shear  simulations  on  polycrystalline  NM. 

•  Investigate  size  dependence  more  thoroughly. 

•  Investigate  time  scale  for  reactions  to  occur  by 
calculating  activation  energies  for  proposed  reactions 
in  NM. 

•  Investigate  hydrostatic  pressure  dependence  of 
results. 
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■  Uncovered  short-time  scale  shear  dynamics  and  possibly 
initial  stages  of  reaction  pathway. 

■  NM  transforms  into  a  disordered  material  under  shear 
stress  at  20  Gpa  and  298  K,  but  no  reactions  are 
observed. 

■  Significant  orientation  dependence  on  the  shear 
response  of  NM:  max  shear  stress  5. 2-1. 6  Gpa. 

■  Shear  stress  characterized  by  the  ability  or  lack  of 
ability  for  NM  molecules  to  rotate  or  planes  of  NM 
molecules  to  shift  with  the  shear  stress. 

■  Response  of  NM  with  a  GB  under  shear  characterized  by 
rotation  of  entire  grain  or  rotation  of  molecules  in  the 
grain. 
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Maximum  Shear  Stresses  /iftL 

Orientation 

<hkl> 

Max  ISxy 

|  (Gpa)  Max  IS^I  (Gpa) 

Max  |Syz|  (Gpa) 

<i-i-i> 

4.21 

5.24 

4.03 

<10-1> 

2.15 

2.70 

4.67 

<11-1> 

3.73 

3.34 

4.54 

<110> 

3.03 

4.53 

2.98 

<010> 

3.92 

3.46 

1.95 

<1-11> 

4.06 

3.03 

2.96 

<111> 

4.08 

3.27 

3.01 

<1-10> 

3.43 

3.46 

2.85 

<01-1> 

4.38 
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